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Based on the density functional theory with hybrid functional approach, we have studied the 
structural and thermodynamic stabilities of Cu2MSnX4 (M = Zn, Mg, and Ca; X = S and Se) 
alloy, and have further investigated the electronic and optical properties of stable Cu2MgSnS4 and 
Cu2MgSnSe4 phases. Thermal stability analysis indicates that Cu2MgSnS4 and Cu2MgSnSe4 are 
thermodynamically stable, while Cu2CaSnS4 and Cu2CaSnSe4 are unstable. The ground state 
configuration of the compound changes from kesterite into stannite structure when Zn atoms are 
substitued by larger Mg or Ca atoms. An energy separation between stannite and kesterite phase 
similar to that of CZTS is observed. Calculated electronic structures and optical properties suggest 
that Cu2MgSnS4 and Cu2MgSnSe4 can be efficient photovoltaic materials. 


I. INTRODUCTION 

Quaternary Cu 2 ZnSnS 4 (CZTS) and Cu 2 ZnSnSe 4 (CZTSe) compounds have drawn significant interests in the past 
years because of the abundance of their constitutiim elements, the environment-friendliness, and the appropriate band 
gaps (1.4 — 1.6 eV and 0.95 — 1.05 eV for CZT^l® and CZTS^^®). They are regarded as a potential alternative 
to Cu(In,Ga)Se 2 (GIGSe), which is the most efficient thin-film solar cell material with a conversion efficiency record 
of 21.7%]^ Mi xing S and Se in Gu 2 ZnSnS 2 ,Se 4 _a; (GZTSSe) can tune the band gap to achieve the best efficiency 
of 12.6% so far.IllJ However, this efficiency is still lower than the requirement of commercially viable device. The 
physical origin of low efficiencies is complex. Gokmen et al. pointed out that the existence of electrostatic potential 
fluctuation in the CZTS(Se) absorber layer could be an important factor that deteriorates the solar cell performances 
and the [Guzn+Zncu] defect cluster could be the origin for the electrostatic potential fluctuation.ll^ One trial to reduce 
the density of [Guzn+Zncu] defect cluster as proposed by Scragg et al. is to design a thermal process with a low 
cooling rate after the high-temperature selenization or sulfurization process to obtain ordered GZTS.I^ To eliminate 
the antisite defects originating from intermixing of Cu and Zn atoms which are actually with very similar atomic 
sizes, an alternative approach is to substitute other elements for Cu or Zn to possibly suppress the defects formation 
and thus improve the sample quality. In particular, the substitution large-size atoms for Cu or Zn atom is suggested 
in these Cu-poor CZTS and CZTSe. 

For the substitution for Zn atom, Cd, Hg, Mg, Ca, and Sr as selected metals have been studied.li^lti^ These designs 
obviously increase the band gap of corresponding sulfides or selenides. From the stability analysis based on the 
generalized gradient approximation (GCA) calculation, the substitutions Cd and Hg for Zn are thermodynamic stable, 
while other substitutions are unstable against the phase-separation into the competing binary and ternary compounds. 
However, Cu 2 MgGeS 4 nanoparticles have recently been synthesized, which shows X-ray diffraction pattern in good 
agreement with that of the theoretically proposed structure.^i^ Particularly, the synthesis of Cu 2 MgSnS 4 naiioparticles 
without coexisting with binary compounds further implies the feasibility of Mg substitution for Zn element For the 
contradiction, checking Wang’s calculations, we conclude that GCA has difficulty in accurately describing the energy 
such as the band gap and Wang et did not include all the related phases to analyze the thermodynamic stability. 
Additionally, both Cd and Hg are not environment-friendly. We, therefore in this work, focusing on large size ions 
Mg and Ca and using more precise method, are to investigate the induced effects by the substitution of Zn to explore 
the feasibility of Cu 2 MSnX 4 (M = Mg and Ca; X = S and Se) as photovoltaic materials. 
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II. COMPUTATIONAL DETAILS 

The calculation is carried out within the framework of the density functional theory (DFT). The Vienna ab initio 
simulation package (VASP^ is employed to optimize the crystal structures and calculate the stability and electronic 
and optical properties. The inner electrons are treated by the projector augmented wave method, while the valence 
electrons are expanded in plane waves with a cutoff energy of 400 eV. For the optimization, a conjugate-gradient 
algorithm is used to relax the ions into their instantaneous ground state. In all calculations, the Monkhorst-Pack 
fc-point grids are generated according to the specified fc-point separation 0.02 and the convergence thresholds are 
set as 10“® eV in energy and 0.005 eV/A in force. 

In the standard DFT, the GGA of Perdew-Burke-Ernzerhof (PBE) versiorP^ is adopted to describe the electronic 
exchange-correlation (XC) interactions. However, the standard DFT usually leads to erroneous descriptions for some 
real systems such as underestimating the band gap. One way of overcoming this deficiency is to use the Heyd-Scuseria- 
Ernzerhof (HSE) hybrid functional, where a part of the nonlocal Hartree-Fock (HF) type exchange is admixed with 
a semilocal XC functional, to give the following expressiorPil 

+ (1 - 

+ ( 1 ) 

where a is the mixing coefficient and v is the screening parameter that controls the decomposition of the Coulomb 
kernel into short-range (SR) and long-range (LR) exchange contributions. In this calculation, the HSE exchange- 
correlation functional is set as 25% mixing of screened HF exchange to PBE functional, namely, a = 0.25. Noticeably, 
all results presented in this paper are obtained from HSE calculations. 


III. RESULTS AND DISCUSSION 

Quaternary I-H-IV-VI semiconductors derived from cation mutation of CuInSe 2 -like I-HI-VI compound may adopt 
one of the three structural phases (Fig. 1), i.e. kesterite structure (KS) with the space group of /4, stannite structure 
(ST) with the space group of /42m, and the primitive mixed CuAu-like structure (PMCA) with the space group 
of PA2m. While the first structure is derived from the chalcopyrite structure, the latter two are derived from the 
CuAu-like structure. The aforementioned structures were considered for the analysis of ground state structure of 
Cu 2 MgSnS 4 (CMTS), Cu 2 MgSnSe 4 (GMTSe), Cu 2 CaSnS 4 (CCTS), and Cu 2 GaSnSe 4 (GGTSe). From calculated 
total energies based on HSE functional (Table I), it is found that PMGA is the least stable structure among these 
three. While we have confirmed the previous calculation results showing that KS is the ground state structure for 
Zn-based compounds.^ On contrary, both Mg-based and Ca-based compounds are stabilized at the ST phase, which 
is similar to other substitution induced by large size atoms such as Cu 2 CdSnA 4 and Cu 2 HgSnA 4 (A = S and Se). 
The energy difference between KS and ST is 47.0 meV per unit cell for GZTS and 80.4 meV per unit cell for CZTSe, 
which are according with Ghen’s results.^ In CMTS and GMTSe, this energy difference between KS and ST is 40.2 
and 42.1 meV per unit cell, respectively. But this energy difference has a large of increase in CCTS and GGTSe. This 
value respectively reaches to 543.8 and 536.0 meV per unit cell for CCTS and GGTSe. The energy difference between 
ST and KS phase of CMTS and GMTSe is comparable to that of GZTS, hence CMTS and GMTSe based on ST phase 
are potentially stable at thermodynamic equilibrium. 

As visualized from the crystal structure, the small total energy difference between the two structures could result 
in the existence of large amount of antisite defects. In Zn-based materials, based on the total energy difference value 
displayed in Table I, at high-temperature, the small total energy difference between the KS and PMCA phases in 
GZTS implies the easy formation of PMCA secondary phase or [Znsn+Snzn] related defect complex. In the Mg-based 


TABLE I: Calculated total energies (meV per unit cell) relative to the most stable phase among KS, ST, and PMCA for 
Cu 2 MSnA 4 (M = Zn, Mg, and Ca; A = S and Se). 


system 

Eks 

Est 

EpMCA 

CZTS 

0 

47.0 

66.4 

CMTS 

40.2 

0 

95.9 

CCTS 

543.8 

0 

694.9 

CZTSe 

0 

80.4 

131.1 

CMTSe 

42.1 

0 

148.5 

CCTSe 

536.0 

0 

754.5 







3 


TABLE II: The optimized lattice constant, the anion displacement parameter u, the nearest atomic distance, and the calculated 
band gap for Cu 2 MSnA 4 (M = Zn, Mg, and Ca; X = S and Se) 


System 

a (A) 

c(A) 

u 

dcu-x (A) 

dM-X (A) 

dsn-X (A) 

E, 

Expt. 

(eV) 

GCA 

HSE 

CZTS (KS) 

5.475 

10.943 

0.25±0.0157 

2.321 (2.325) 

2.370 

2.471 

1.4-1.6“ 

0.09 

1.60 

CMTS (ST) 

5.568 

10.932 

0.25±0.0006 

2.330 

2.471 

2.464 


0.19 

1.80 

CCTS (ST) 

5.903 

10.483 

0.25±0.0206 

2.363 

2.726 

2.449 

— 

0.55 

2.13 

CZTSe (KS) 

5.771 

11.537 

0.25±0.0183 

2.437 (2.444) 

2.499 

2.623 

0.95-1.05“ 

0.02 

1.00 

CMTSe (ST) 

5.864 

11.562 

0.25±0.0006 

2.447 

2.611 

2.617 

— 

-0.05 

1.16 

CCTSe (ST) 

6.173 

11.186 

0.25±0.0185 

2.476 

2.859 

2.602 

— 

0.18 

1.62 


“Ref. [THi 
'’Ref. [T8| 
“Ref. CHI 


TABLE III: Listing of copper chemical potential limit and volume of stable chemical potential range for Cu 2 MSnS 4 and 
Cu 2 MSnSe 4 (M = Zn, Mg, and Ca) showing that replacement of Zn by Mg forms energetically stable compound. 



M = Zn 

M = Mg 

M = Ca 

Quaternary sulfide /rcu range (eV) 

-0.72/0 

-0.67/0 

Unstable 

Quaternary selenide /icu range (eV) 

-0.65/0 

-0.62/-0.01 

Unstable 

Volume of stable chemical potential of Cu 2 MSnS 4 

0.076 

0.038 

— 

Volume of stable chemical potential of Cu 2 MSnSe 4 

0.057 

0.023 

— 


materials, the PMCA secondary phases can not easily form both in CMTS and CMTSe because of the relatively large 
total energy difference between the ST and PMCA structures. In Ca-based materials, the much larger atomic radius 
of Ca than that of Cu can induce a strong mechanical stress and can suppress the [Cuca+Cacu] complexes. The total 
energy difference between ST and PMCA also increases so that the PMCA secondary phases and [Casn+Snca] defect 
complexes can be suppressed. 

For the ground-state structure of Cu 2 MSnA 4 (M = Zn, Mg, and Ca; A = S and Se), the lattice constants are 
listed in Table II to compare the substitution effects. The substitution of Mg for Zn leads to expansions along a 
and b directions. From the atomic distances, the bond lengths of dcu-x and d-M-x is found to increase with M 
atomic radius. Except for CMTSe, the substitution induces a shrink along c direction, which is possibly related to the 
shortening of bond length of dsn-x in Mg- and Ca-based compounds. Noticeably, the anion (S or Se) displacement 
parameter u deviates from the ideal value of 0.25 for Zn- and Ca-based materials, while it is almost ideal (u = 0.25) 
for Mg-based compounds. 

At experiment, more stable secondary and ternary impurity phases often decrease the crystallizing quality of 
quaternary compounds, or even no stable quaternary phase is obtained. So we have also examined the thermodynamic 
stability of the materials above in order to provide a guidance for their synthesis. The requirement that the Cu 2 MSnA 4 
(M = Zn, Mg, and Ca; A = S and Se) phase is more stable than any secondary phases at thermodynamic equilibrium 
imposes a restriction on the range of chemical potentials. Satisfying the constraint requires that 

AiJf(Cu2MSnA4) 

= 2A/rcu + + A^Sn + 4A^x, (2) 

and for any secondary phases 

Ai?f(secondary phases) > (3) 

i 

where AH{ denotes the formation enthalpy of the corresponding phase, rii and Afj^i denote the number of atoms and 
the chemical potential of the corresponding element Constrained by Eq. (2), only 3 out of 4 chemical potentials 
are independent. We hence choose by convention the chemical potentials of metallic elements to be independent 
variables. It is clear that any Afj, should have a value less than 0 to avoid segregation of elements. Binary and ternary 
compounds further restrict the range of stable chemical potential. In the calculation, element Cu, M =Zn/Mg/Ca, Sn 
and A =S/Se, binary metal sulfides and selenides CU 2 A, CuA, MX, SnA, SnA 2 , and ternary compounds Cu 2 SnA^^, 
Cu 3 SnA 4 l^are the secondary phases that have been considered. 

We adopt the range of copper chemical potential and the volume of stable chemical potential as an indicator of 
the feasibility of growing the Cu 2 MSnA 4 material as copper vacancy is likely a candidate of p-type conductivity. 
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The ternary phase Cu 3 SnX 4 is not involved in the boundary of stable chemical potential range when HSE exchange- 
correlation function is used. The copper chemical potential limits and volume of allowed chemical potential range are 
shown in Table III. The copper limit of Mg-based material is slightly smaller than that of Zn-based counterpart, with 
the volume of allowed chemical potential remaining 40—50% of Zn-based compound, and Ca-based material is unstable 
with respect to its secondary phases. The mixing of Cu 2 SnX 3 and MgX to form Cu 2 MgSnX 4 is thermodynamically 
favorable if 


Ai/f(Cu2MgSnX4) < Ai7f(MgA) -t Ai7f(Cu2SnX3) (4) 

even with entropy effect not considered. The reason that previous worlP^ using GGA functional had not reveal 
the thermodynamic stable range of Gu 2 MgSnA 4 is mainly due to GGAs inability to correctly evaluate ground state 
energy of Gu 2 MgSnS 4 and Gu 2 SnS 3 . For example, a comparison of GGA result to that of HSE shows 0.98 eV rise in 
formation energy in Cu 2 MgSnS 4 , while the rise of formation energy for Gu 2 SnS 3 and MgS is 0.62 eV and 0.25 eV. This 
subtle energy difference actually resulted in a significant change in the phase diagram as it completely eliminates the 
stable chemical potential region. The relative smaller volume of allowed chemical potential for Mg-based compound 
is likely due to the relatively larger internal stress of the material. Since ZnS (ZnSe), Gu 2 SnS 3 (Gu 2 SnSe 3 ) and GZTS 
(CZTSe) have similar lattice constants, the bond length between Zn and S (Se) in CZTS(Se) is almost identical to that 
in the binary ZnS (ZnSe) compounds. As a result, internal stress around Zn elements is expected to be small in these 
compounds. In contrast, the ionicity of Mg and Ga results in an energetically more stable halite-structured MgS(Se) 
and GaS(Se), with significantly longer bond length and higher coordination number than the zinc-blende counterpart. 
Atoms in Mg-based compounds are likely to be more stressed. This internal stress may lead to a relatively narrower 
stable chemical potential range for GMTS and CMTSe. 

To evaluate the potential of CMTS and GMTSe as photovoltaic materials, however, we have further investigate 
their electronic and optical properties and compare them with that of Zn-based materials. As shown in Table II, 
GGA strongly underestimates the band gap. Adopting HSE hybrid functional, we reproduce the experimental band 
gaps of GZTS and GZTS^i®. ST structured GMTS and GMTSe are direct band gap materials, with calculated band 
gap as 1.80 eV and 1.16 eV respectively, which are correspondingly larger than those of GZTS and GZTSe. Ga-based 
materials have bigger band gap of 2.13 eV for GGTS and 1.62 eV for GGTSe, which mainly comes from the stronger 
hybridization interaction between Sn-s and X-p states as well as from the weaker coupling interaction between M — d 
and X — p states. Seen from the band structures and the projected density of states (PDOS) shown in Fig. 3, Mg 
substitution does not result in a change in nature of valence and conduction band and makes no contribution to the 
electronic states near the Fermi level. The valence bands mainly come from the hybridization of Cu—3d and S—3p 
(Se—4p) electronic states and the first-group conduction bands mainly come from the mixture of Sn—5s and S—3p 
(Se—4p) electronic states. The small increase of band gap in Mg-based materials mainly results from the enhancement 
of interaction between Sn and S (Se) induced by the bond length shortening comparing to Zn-based materials. 

The optical properties of GMTS and CMTSe are compared with those of Zn-based materials. Here, the calculation 
of the optical properties is based on the complex dielectric function s{uj) = ei(uj) +ie 2 {uj) describing the polarization 
response of the material to an externally applied electric fielcP^. The dielectric function of CMTS and CMTSe are 
shown in Fig. 4 and correspondingly compared with GZTS and CZTSe. The onset to the response of the imaginary 
part e 2 {oj) of the dielectric function reflects the band-gap energy, and in the vicinity of this onset of the ST CMTS 
and CMTSe show strong anisotropy compared to the KS GZTS and CZTSe. This is a consequence of the 

stronger and positive crystal-field splitting in the ST compounds. At somewhat larger energies (in the regions about 
2.0 — 4.0 eV for sulfides, 1.5 — 3.5 eV for selenides) opposite anisotropy appears: £f^(w) < er(w) for ST CMTS and 
CMTSe, and £f’^(a;) > £™(a;) for KS GZTS and CZTSe. The anisotropy is reflected also in the real part £i(a;) of the 
dielectric function. The selenides show larger high-frequency dielectric constant £oo in accordance with the smaller 
band gap Eg. The ST Mg-based compounds show slightly smaller high-frequency dielectric constant than the KS 
Zn-based materials. The value of £oo is estimated to be about 3.9 and 4.8 for CMTS and CMTSe, respectively. 

The absorption coefficient a{u}) for ST CMTS and CMTSe, together with those for GZTS and CZTSe, for the 
optical absorption from the valence bands to the lowest conduction bands are shown in Fig. 5. Since the absorption 
is obtained directly from the dielectric functiorPSI, similarities in the polarization response are reflected also in the 
absorption coefficient. Thus, all four compounds have comparable absorption, although with different photon energy 
for the onset. A large band-edge absorption coefficient of ~ 10® cm“^ is observed for all the four materials. 


IV. CONCLUSIONS 

Using the first-principles calculations and combining with the HSE hybrid functional, we have studied the induced 
effects by the substitution of Zn in Cu 2 ZnSnA 4 {X = S and Se), from structural stability, thermodynamic stability, 
electronic characteristics, and optical properties. Large size atom Mg and Ca substituting for Zn, we found that 


5 


the ground-state structure is stabilized at ST configuration instead of KS as Zn-based compounds. The result of 
total energy difference between KS and ST or between ST and PMCA indicates that the substitution larger size Ca 
atom for Zn atom can efficiently suppress the antisite defects. Unfortunately, Ca-based materials are not thermody¬ 
namically stable with respect to the secondary phases. Mg-based materials are thermodynamic stable, with energy 
difference between the ST and KS phases comparable to that of CZTS and CZTSe. Additionally, [Mgsnd-SnMg] defect 
complexes are expected to be largely reduced in CMTS and CMTSe. Mg substitution does not change the electronic 
characteristics but opens a slightly bigger band gap than Zn-based compounds. The optical absorption spectra show 
the feasibility of Cu 2 MgSnS 4 and Cu 2 MgSnSe 4 as a photovoltaic material. 
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0Cu Ozn/Mg/Ca O S/Se 



FIG. 1: (Color online) Three possible configurations of Cu 2 MSnX 4 (M = Zn, Mg, and Ca; X = S and Se): KS, ST, and 
PMC A. The purple arrows mark the [Cuzn+Zncu] defects complexes that result in the mutual transformation between KS and 
ST structures. The black arrows mark the [Znsn+Snzn] defects complexes that result in the mutual transformation between 
ST and PMCA structures. 
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FIG. 2: (Color online) Illustration of phase diagram of CMTS and CMTSe. The large surface bound by Cu 2 SnS 3 (Cu 2 SnSe 3 ) 
is the front surface and that of MgS (MgSe) is the back surface. Indicators of Cu-chemical potential are plotted at Cu-richest 
limit (A-D), Cu-mild range (E-F) and Cu-poorest limit (G). 
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FIG. 3: (Color online) Calculated band structures and projected density of states (PDOS) for ST Cu 2 MgSnS 4 (top) and 
Cu2MgSnSe4 (bottom). 
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FIG. 4: (Color online) The dielectric function e{uj) = ei(tu) + ie 2 {oj) of Cu 2 MgSnS 4 and Cu 2 MgSnSe 4 . The red line represents 
the real part ei(i^) and the blue line represents the imaginary part The dielectric function is divided into the transverse 

contribution (xx; solid lines) and longitudinal contribution (zz; dotted lines). For comparison, the dielectric function of KS 
Cu 2 ZnSnS 4 and Cu 2 ZnSnSe 4 are also presented. 
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FIG. 5: (Color online) Calculated absorption coefficient a{u}) in visible light region of ST Cu 2 MgSnS 4 (a) and Cu 2 MgSnSe 4 
(b). For comparison, the absorption coefficient a{u}) of KS Cu 2 ZnSnS 4 and Cu 2 ZnSnSe 4 is also presented. 








